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Abstract 

In the framework of quantum gravity propagating in large extra dimensions, we an- 
alyze the inclusive radiative emission of Kaluza-Klein spin-2 gravitons in the two- 
fermions decays of massive gauge bosons, heavy quarks, Higgs bosons, and in the 
two-massive gauge bosons decay of Higgs bosons. We provide analytical expressions 
for the square modulus of amplitudes summed over polarizations, and numerical re- 
sults for the widths and branching ratios. The corresponding decays in the Z, top 
quark, and Higgs boson sectors of the standard model are analyzed in the light of 
present and future experiments. 



1 Introduction 

After the recent proposal of Arkani-Hamed, Dimopoulos, and Dvali (ADD) on quan- 
tum gravity propagating in large extra dimensions JT], there has been an intense 
theoretical activity on this subject [Qj-0. In [I], it was pointed out that if compact- 
ified extra dimensions exist, with only gravity propagating in the bulk and standard 
matter with gauge fields confined in the usual 3+1 dimensional space, then the fun- 
damental scale of quantum gravity could be much lower than the Planck scale Mp. 
In particular the weakness of gravity might be due to the large size of the compact- 
ified extra dimensional space.0 Indeed, in this scenario, the Newton constant Gat in 
* For a realization of large extra dimension scenarios in the framework of string theories, see j?]]. 
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the 3+1 dimensional space is related to the corresponding Planck scale Mr, in the 
D = 4 + 5 dimensional space, by 



G N X = 8tiR 5 M 2 d +5 (1) 

where R is the radius of the compact manifold assumed here to be on a torus. 

Large extra dimensions can therefore provide a new solution to the hierarchy prob- 
lem and open new attractive scenarios |lj . In particular, if Mr, ~ TeV then deviations 
from the Newton law are expected at distances of order R < 10 32//<5_19 meters ||. The 
present experimental sensitivity in gravity tests is above the millimeter scale, and the 
solution to Eq.(U), with M D ~ TeV, requires 5 > 2. A dramatic phenomenological 
consequence of this theory is that quantum gravity effects could be sizeable already 
at the TeV scale, and could be tested at present and future collider experiments. 

After integrating out compact extra dimensions, the Einstein equations in the four- 
dimensional space describe massive Kaluza-Klein (KK) excitations of the standard 
graviton field.[| These KK excitations are very narrowly spaced in comparison to the 
M D scale, the mass splitting (Am G ) being of order Am G ~ = Mr fff^) , where 
the reduced Plank mass M P is defined as Mj, = (%kG n )- 1 . The couplings of the KK 
gravitons to the standard matter and gauge fields is therefore universal and equal to 
their zero modes, and hence suppressed by l/Mp. On the other hand, in the case 
of inclusive production (or virtual exchange) of KK gravitons, remarkably, the sum 
over the allowed tower of KK states (which could be approximated by a continuos) 
gives a very large number. This number exactly cancels the suppression factor 

( F \ 2+S F 

associated to a single graviton production, replacing it by [jf~) , where E is the 
typical energy of the process. Therefore, if Mr, is in the TeV range, quantum gravity 
effects might become accessible at future collider experiments. 

Another interesting possibility for the solution of the hierarchy problem, as sug- 
gested by Randall and Sundrum ||, is to have a non-factorizable geometry where the 
4-dimensional massless graviton field is localized away from the brane where standard 
matter and gauge fields live. The main signature of this scenario is quite different 



from the one arising from the ADD scenario in collider experiments ||10|| . Indeed, 
widely separated and narrow spin-2 graviton modes are expected. 

In the present paper, we restrict our analysis to the ADD scenario. In this frame- 
work, the relevant physical processes in e + e~ and hadron collider experiments have 
been first analyzed in |JJ (see also [f|,f|). They can be classified in: a) direct pro- 
duction of KK gravitons and b) virtual gravitons exchange. In the first case, the 
best signatures corresponding to the final state would be a photon associated with 



t For a detailed discussion about the effective four-dimensional theory, see . 
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missing energy (in electron colliders) or jet + missing energy (in hadron colliders). 
In the latter case, the gravitons exchange will induce local dimension-eight operator 
(associated with the square of the energy momentum tensor) that will affect the stan- 
dard four fermion interactions processes. The main conclusion of || is that searches 
at LEP2 and Tevatron can probe the fundamental scale up to approximately 1 
TeV, while the CERN Large Hadron Collider (LHC) and linear e + e~ colliders will be 
able to perturbatively probe this scale up to several TeV's. 

In the present scenario, for any new heavy particle with mass close to Md ~TeV, 
the gravitational radiation induced in its decays might become important .[] Indeed, 
the suppression factor in the branching ratio will be given in this case by ( J , 
where M is the mass of the decaying state |IJ. For instance, new particles at the TeV 
scale are expected in some models where the grand-unification scale is lowered down 
to the TeV scale by the appearance of new compact extra dimensions where standard 
model (SM) fields live [|T^]-|]TB|. Such extra dimensions are a natural consequence 



of string theories with large radius compactification. These scenarios could provide 
both a natural explanation for the fermion mass hierarchy (since the fermion masses 
evolve with the mass scale by a power law dependence |I3[), and a natural higher- 



dimensional seesaw mechanism for giving masses to light neutrinos Jl5|. Therefore, 
in a unified picture of gauge and gravitational interactions with unification occurring 
at around the TeV scale, we should expect new Kaluza-Klein excitations of the SM 
particles at the TeV scale. In the decays of such states, the gravitational radiation 
could give rise to relevant effects. 

The aim of the present work is to provide, in the framework of gravity propagating 
in large extra dimensions, analytical and numerical results for both differential widths 
and inclusive branching ratios of gravitational decays of heavy particles, versus the 
number of extra dimensions. In particular, we consider the following classes of decays 

V. // - /, /,-.,- • G. f^f^V + G, H^VV + G (2) 

where V, H, and represent a generic massive gauge boson, Higgs boson, and fermion 
field (with i ^ j), respectively. We retain the masses of all the particles in the final 
states, except in the decay —>■ fj V + G where the final fermion is assumed massless. 
We then apply our results to the analysis of the Z, W, top quark, and Higgs boson 
decays in the SM. On the other hand, our results can be easily applied to more general 
cases, too. 

We restrict our discussion to the spin-2 gravitons, and do not include the cor- 
responding decay modes into scalar gravitons (graviscalars, with J=0) since their 

tThis does not include scenarios with brane deformations, see for instance Refs.|ll| where 
KK(n) tower states of SM particles can have tree-level 2-body decays in KK{n) — > KK(n — 1) + G. 
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amplitudes get smaller with respect to the J=2 ones, being suppressed by a term 



proportional to u = 1/^/3(5 + 2)/2 0. A special case is provided by scenarios where 
the Higgs boson can have a mixing to graviscalar field through the coupling to the 
Ricci scalar f|,||. In these scenarios, the inclusive decay of the Higgs boson in all 
the allowed tower of KK graviscalars is very large, and leads in practice to a sizeable 
invisible width for the Higgs boson ||. 

In the framework of gravity propagating in large extra dimensions, in [|17[] the 
decay Z — > / / + G has been analyzed for both J=2 and J=0 gravitons, versus high 
precision LEP1 Z-pole data. Only numerical results are provided for J=2. As shown 
in section 3, our results for the inclusive total width T(Z — > f f + G) agree with JL7|. 

The paper is organized as follows. In section 2, we define the interacting lagrangian 
describing massive gauge bosons coupled to fermions and Higgs fields, and the cor- 
responding energy momentum tensor which enters the coupling to the graviton field. 
In section 3,4, and 5, we give the analytical and numerical results for widths and 
branching ratios, and discuss the corresponding decays for the SM Z/W, top quark, 
and Higgs boson, respectively. In section 6, we present our conclusions. In appendix 
Al, we report the relevant Feynman rules for the gravitational interaction vertices, 
and in appendix A2 we give the analytical expressions for the square modulus of the 
amplitudes. 



2 Effective Lagrangian 

The coupling of gravity with standard matter and gauge fields in D-dimensional space 
is given by the lagrangian Co @ 

^D = —^HT AB h AB , A =(//,«), /i = 0,...,3, i = 4,...,D-l (3) 

where M^ +s = Mj,(2irR)- 5 , T AB is the energy momentum tensor, h AB the graviton 
field in a D-dimensional space, and the A and B indices refer to the D-dimensional 
space. The sector of the energy-momentum tensor Tab containing standard matter 
and gauge fields is assumed here to have non-zero component only along the A, B = 
H, v directions.^ After integrating out the compactified extra dimensions in the D- 
dimensional action, the resulting (effective) four dimensional theory is described by 
KK graviton fields h> n ' ^ which have the same universal coupling to the SM particles 
as their massless zero-mode (n=0) 0. Then, in four dimensional space, the effective 

§ This can be realized assuming that SM particles correspond to brane excitations and the brane 
itself does not oscillate in the extra dimensions. 
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Lagrangian is given by 

£e// = ^£T^)^, (4) 

where Mp is the reduced Planck mass, and T^ v is the SM energy momentum tensor. 

In this section, we fix our conventions for the Lagrangian £ and its energy momen- 
tum tensor T^ v which are relevant for the processes we are considering. In particular, 
we generalize the fermion fields (/j) couplings of the SM in the weak gauge boson 
(V) and Higgs (H) sectors. This parametrization might be particularly useful in a 
generalization of the SM interactions including KK excitations of SM fields, when SM 
fields are assumed to propagate in other extra-dimensions. 

In Minkowski space, after spontaneous gauge symmetry breaking, the relevant 
Lagrangian in the unitary gauge is given by 

£ = £ F + £ v + £ H + h.c. 

1 M 2 
A/ = — A F, V {V)F^{V) + -2-V ll V 

C H = \ (d,H d»H - M 2 H H 2 ) + J2 \ (A fi H) + 9 -M v (V^ H) (5) 

where 

= l -d»^-ig{g v + g Alh )K^V, 
F^(V) = d^Vy — d v Vn (6) 

dfj. = — dfi, and gy t A represent axial and vectorial couplings (e.g., in the W case 
gy = —gA = ^75)- Kij is a unitary matrix which in this case generalizes the usual 
CKM matrix. Notice that we have restricted our Lagrangian to describe only abelian 
gauge bosons, since we will consider processes involving at most two gauge bosons in 
each interaction. The standard W and Z couplings to the Higgs boson and fermions 
can be easily recovered by this lagrangian. 

In order to obtain the expression for the energy momentum tensor in eq.(^), 
it is useful to rewrite £ in general space-time coordinates with the metric g^". As 
usual, when there are fermion fields, this is simply achieved by the following standard 
procedure. The Minkowski metric if" is replaced by the general metric if — > g^ v 
expressed in terms of the Vierbein fields ej 1 (i.e., g^ u = Y,a e ^ e a 1 where a and /1, v 
are the Minkowski and world indices, respectively) inside Eq.(|5]), and the lagrangian 
£ is multiplied by \f—g, where g is the determinant of g^ u . Then, the expression for 
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Tfa, can be derived by expanding ej 1 around the flat metric 

At the first order in the h^ v expansion, T^ v is given by 

T, v = I/, (j^j + 7 „Dj?) fj - r,^ (irfl% - Mi5 ij ) fj 
+ F m F\ + MyV^Vy + ^ F aP - ^M v V a V°^ 

+ d^H d v H + gM v V[j~V v H - -rj^ (d a Hd a H - M 2 H ) 

+ V^(xJifiH + 9 -M v V a V a H} (8) 

where the sum over the i,j fermion flavours is assumed. Notice that, at first order 
in hj 1 , there is no distinction between latin (a) and greek indices (//), being all the 
contractions performed by rj^, and g^ u = rf + ~- (W + h ufl ) + 0{h 2 ). 

By inserting eq.(§) in eq.(^), the corresponding Feynman rules for interaction 
vertices with both three-line and four-line vertices (and only one graviton emission) 
are easily obtained. We report their expressions in Appendix Al. All the four-line 
vertices contribute to the matrix elements involving on-shell spin-2 fields except in 
the case of the fermion-Higgs couplings. In the latter case, the HffG vertex is 
proportional to the trace of the energy momentum tensor (see eq.([S])) and, therefore, 
its spin-2 component vanishes. 



3 Heavy Gauge bosons and Z /W decays 

We start our analysis by considering the following decay 

V( Pv )-Hpi)f l AP2)+G( PG ) (9) 

where V in this case can be both a U(l) massive gauge boson and a non-abelian 
SU(N) massive gauge boson, G is a massive spin-2 field, and fi is a generic fermion 
field. We also set = 5ij in eq.(^). The four Feynman diagrams relevant to this 
process are shown in figs. |l|a-|l]d. We see that the diagrams in figs, la-lc are obtained 
by attaching the spin-2 field in all possible ways to the external legs of the main 
diagram V —> f f, while fig. Id is given by the contact term VVHG. 

Since we are interested in studying unpolarized processes, we recall here the for- 
mula for the sum over polarizations in the case of a massive spin-2 field. This is given 
by 

5 

^{k, cr) e af3 (k, a) = P^^k) (10) 

(7=1 
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(a) (b) (c) (d) 

Figure 1: Feynman diagrams (a-d) for the decays V — > / / + G or /, — > fjV + G . Analogous 
diagrams in the Higgs sector can be simply obtained by these diagrams, by replacing V —* H 
for the decay H — > / / + G, and also replacing (/, /) — > (V, V) for the decay if — ► T^F + G. 
In the case H —> ff + G the (d) diagram vanishes, as explained in the text. 

- 7T-T (V^ak^kp + r] u(3 k^k a + ri^Kkp + "q^k^kp) 

+ I (vi* + -Kk^kA (rjap + \k a kA , (11) 

6 V m G J \ m G J 

where e Atl/ (fc,a), uig and k are the polarization tensor, mass and momentum of the 
spin-2 field, respectively, and the index a runs over the polarization states. Note 
that the projector P^ ua p, which is symmetric and traceless in both (//, v) and (at, 0) 
indices, satisfies the transversality conditions k^P^ ua p = k a P lxva i3 = 0. Then, by using 
the Lagrangian in eqs.(^) and (|j), the square modulus of the amplitude M. summed 
over all final polarizations and averaged over the initial ones, is 

pol P 

(12) 

where g is the gauge coupling constant, and gy,A are the vectorial and axial couplings 
defined in eq.(|5[). For fermion /, Nf represents the sum over quantum numbers, whose 
generators commute with the gauge group generator associated to the vector V, like, 
for instance, the fermion color number in the case of realistic Z or W decays. We 
assume the two fermion masses degenerate (i.e., Mf i = Mf j = M/). Then, we define 
the Mandelstam variables t and u as 

t = jp (Pi+Pg) 2 ~ Xf, u= j^(p2+Pg) 2 ~ x f , s = x G -t-u (13) 
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where Xf = -^f , xg = and My, ma, Mf are the masses of the gauge boson, 
graviton, fermion, respectively. The analytical expressions for the functions Fy^ (t, u) 
(also depending on the variables Xf and xq) can be found in the appendix A2. 

It is worth noticing that, despite the presence of l/m G terms in the sum over 
polarizations for the massive spin-2 fields, in the final expression for F y ±} (t,u) (and 
analogously for the other decay functions in the appendix A2) mg appears only with 
positive powers. The cancellation of l/m G terms in the total amplitude is indeed 
ensured by the conservation (at the zeroth order in h^ v ) of the on-shell matrix elements 
of the energy momentum tensor in eq.(H). Therefore, the terms proportional to in 
eq.flnp do not contribute to the final amplitude. As a severe check of our results^, we 
used the complete expression for P^ap in eq. ( |TT| ) an d explicitly verified this property. 

Although the limit for mc — ^ is smooth, our results for the square amplitudes 
summed over polarizations is not supposed to coincide in this limit with the massless 
graviton contribution. This is due to the well-known van Dam-Veltman discontinuity 
|T8| . Our results only hold for vtlq ^ 0. Indeed, the emission of a massless graviton 
should be calculated by using the proper massless projector (see, e.g., fl8fl), that differs 
from the massive one.[] For the purpose of our analysis the effect of not taking into 
account this discontinuity is not relevant, since the square modulus of the amplitude 
for a single massless graviton contribution is suppressed by 1/Mj,. 

As said above, we are interested in analyzing inclusive processes, where one sums 
over all the kinematically allowed KK graviton states. The mass splitting Aitlg 
between different excitations is given by 

Am °~R= M °U) (15) 

(e.g., in the case 5 < 4, for Mo ~1 TeV, Arric is less than a few KeV's). This allows 
one to approximate the KK modes as a continuous spectrum, with a number density 

^The functions appearing in appendix A2 were computed by FORM [19] ■ 
II In particular, P^ a (3{k) in eq.(|ll|) for the massless case becomes Jl8| 

P ^ a p{ k ) = \ {V^cxVufi + VuPVva ~ VlfVap) + ■■■ (14) 

where the dots stand for any term containing at least one graviton momentum. Notice that, also 
in this case, the terms proportional to the graviton momentum give zero when contracted with the 
on-shell matrix elements of T M „, due to the conservation of the energy- momentum tensor. Therefore, 
the discontinuity in the limit ttlq — > arises from terms that do not contain the graviton momentum 
in the two projectors. In particular, the only difference in the relevant terms of eqs.(|l4|) and ( |Tl| ) is 
given by the coefficients of ?7 m „77q/3, that, when contracted with the on-shell matrix elements of the 
energy momentum tensor, give terms proportional to the trace Tff. Since we are analyzing massive 
particles, the trace of does not vanish, so in the limit ^Owe should expect a discontinuity. 



of modes (dN) between mg and mg + dm G given by 0] 

M 2 

dN = Ss-i^^dma. (16) 

Here, S^-i is the surface of a unit-radius sphere in 8 dimensions.!**] Then, the inte- 
gration over the number of KK states cancels the factor 1/Mp of the single graviton 
emission. 

Finally, the result for the inclusive total width T(V — > ff + Gx), where Gx 
indicates any KK graviton excitation up to the My scale, is given by 

+ (\9v\ 2 -\9A\ 2 )l l v\x A ,x f ,6)) (17) 



where 



I Y h) (x A ,x f ,S) = / dx G (xg)^ 1 / dt / duFy ±} (t,u) (18) 



and 



(— 2x/ + x G - t) t + x G + t ± A 

M± ~ W+^7) 

A 2 

(x G - t) 2 - Ax f x G ) (1 + t 2 - 2 (2xy + t)), x A = (19) 

My 

We defined GV = g 2 /(A\/2M v ) extending the definition of the standard Fermi con- 
stant Gp. A exp is the experimental resolution on the invariant mass of missing energy. 
Notice that the function Iy(x/±,Xf,8), that is proportional to the fermion masses, 
vanishes in the Xf — > limit. The integral over the phase space is also regular in the 
xa — > limit. This property is connected to the absence of graviton mass singularities 
in the square amplitude. 

In table 1, we report some numerical values of the integrals Iy(x^,,Xf,8) in 
eq. (|T8|) , evaluated at ia = and Xf — 0, for the representative cases 8 = 2,3, ... ,6. 
In fig. H, we plot the integrals Iy(x&, Xf, 8) (evaluated at x<\ = and divided by 
their values at Xf — [see table Q) versus Xf, and for 8 = 2, 3, . . . , 6. In fig. we 
plot the differential widths versus the ratio r G = m G /My, and for 8 = 2, 3, . . . , 6. In 
particular, we plot the distribution R, normalized as 

1 dT 

R =^r~^r (2°) 

1 TTlflY til CI 



dr 



a 



**From @, we get S S -i = 2n n /(n - 1)! and S S -i = 2ir n / UtJ)( k + |) for 5 = 2n&ndS = 2n + 1, 
with n integer, respectively. 
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5 


A+) 
1 V 


h 




T v 

A H 


2 


1.20 x 10- 1 


5.53 x 10- 2 


3.25 x 10" 1 


7.78 x 10~ 2 


3 


2.12 x 10~ 2 


8.72 x 10~ 3 


4.95 x 10" 2 


1.17 x 10" 2 


4 


5.85 x 10~ 3 


2.10 x 10~ 3 


1.14 x 10~ 2 


2.67 x 10~ 3 


5 


2.11 x 10~ 3 


6.57 x 10~ 4 


3.35 x 10~ 3 


7.86 x 10" 4 


6 


9.13 x 10~ 4 


2.45 x 10~ 4 


1.18 x 10~ 3 


2.76 x 10" 4 



Table 1: Numerical values for I v ' (xa, Xf, 5), If(xA,Xy,5), ljj(xA,Xf,5), and 
Ih(xa, xy, 8) evaluated at xa = xy = %f = 0, and for 5 = 2,3, 4, 5, 6. 



where J^| ma x stands for the maximum of versus r G . The shape of this distribution 
provides information on the typical fraction of missing energy, due to the KK gravitons 
emission, expected in the decay. We analyze two representative cases: Mf = and 
Mf = 0.2 My. From these results, we see that the position of the maximum (r max ) of 
the distributions R is quite sensitive to the number of extra dimensions, going from 
r max ~ 0.1, for 5 = 2, up to r max ~ 0.5, for 6 = 6. When the mass for the final 
fermion is taken into account, these curves shift toward lower tq values, due to the 
phase-space reduction of the allowed xq range. 

We recall that our perturbative treatment is bound to be valid for decaying par- 
ticles not heavier than Mq. Indeed, being gravity directly coupled to the energy 
momentum tensor, the validity of the perturbative expansion strongly depends on 
the energy scale of the process with respect to the Planck mass Mj> Since in the rel- 
evant scenarios Mo could be close to the TeV scale, this question is not just academic. 
In [0, when considering direct graviton production at colliders, upper bounds on the 
center of mass energy as a function of Mp and number of extra dimensions have 
been obtained by requiring unitarity of the tree-level cross sections for single graviton 
production. In our case, if new particles exist with masses either larger than or close 
to the fundamental Planck scale Mp, then non-perturbative gravitational phenom- 
ena should sizeably affect their decay width. For instance, mass upper bounds that 
somewhat limit the perturbative regime can be obtained by requiring that the rate 
for one graviton emission from the main particle decay does not exceed the rate of its 
main decay. In particular, one can impose 

r(v - //) * (n> 

for any fermion / and for A exp = 0. In the approximation of massless fermions the 
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width of V — > // is given by 



r<v - //) = ^^f 1 (M 2 + M 2 ) (22) 



then from Eq.(^TJ), one obtains 



My < M D ( TFy - :rT77 — ) • (23) 



64 7T 2 \ 2+i 

jr ; (^A,0, 5) S s - 
From the values of the integrals given in table 1, for xa = we find 



M v <M D x (5.4, 4.7, 4.2, 3.8, 3.5) , (24) 

for 5 = 2,3,4,5,6, respectively. For massive fermions, the value of /j> + ^ would be 
smaller (due to a smaller available phase space), giving less stringent upper bounds 
on M v . 

As a phenomenological application of this study, we analyze now the constraints 
on Md which come from negative searches of extra missing energy events in the Z 
decays. In particular, we will consider the process 

Z // + G x (25) 

where / may indicate leptons and quarks. The massless limit for the final fermion 
states is quite accurate in this case. From eq. (|22|) , the decay width of Z — > fafe for 
massless fermions is simply 

V(Z^l}) = N, 2 -f^{\ 9 ' v ? + \ 9 ' A ?) (26) 

where Nf = 1 and 3 for leptons and quarks respectively, and g* A = T% /2, g v = 
{rl — 2Qf sin 6w 2 ^j /2, being T% and the eigenvalues of the third component of 
isotopic spin (T 3 e = — |) and electric charge, respectively. Here, the vectorial and 
axial couplings have been properly normalized after the introduction of the true Fermi 
constant Gp. Then, the branching ratio (BR) for the inclusive decay Z — ► J2f f f + 
Gx, where / stands for any quark or lepton in the final state, is given by 



2+8 



BR(Z - £ // + G X ) = ^± J«(* Aj 0, 5) (27) 

In the case 5 = 2, and for ia = = 0, we obtain 



m z^ff + G x ) = %^ (2. 
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The result in eg. (^8|) agrees with the corresponding one in |L7|, after identifying 
= 2Mf) for 5 = 2, being the definition of M* in [T7J different from Mp in eq.([I|) 
(G" 1 = 4nR s M? +s ). Note that in |7 some eraviscalar contribution is included. 

These results can be extended to the decay W ± — > ff + G, in the massless limit 
for final fermions. In the latter case, the width and inclusive BR can be simply 
obtained from the Z case, by replacing Mz — > My/ in eq. (|27|) . In the case 6 = 2, and 
for xa = 0, we obtain 

5.0 x 1(T 8 



BR(W Z 



/,/' 



M^(TeV) 



(29) 



where f ^ f run over the leptons and quarks with T% = 1/2 and T% = — 1/2. 

The result in eq. (|28|) can be applied to the LEP1 data on Z — > // + E miss , 
corresponding to about 2 x 10 7 Z decays. The SM background, given by the four- 
fermion decay Z — > ffvv, is very small. A few events that are quite in agreement 
with the SM prediction were observed p0| . Assuming, one can push the limit on 
unexpected signals down to BR(Z — ► Z)/ // + Gx) < 10 -7 , one then gets, for 5 = 2, 

M D > 951 GeV (30) 

This limit is not far from what is obtained from the negative searches at LEP2 in 
the channel e + e~ — > 7 + E miss , (where E miss is the missing energy due to gravitons 



emission) and from virtual gravitons effects [21 



4 Heavy Fermions and Top decays 

We consider now a heavy fermion fa decaying into a lighter fermion fj plus a massive 
vector boson and a graviton 

fi(pi)^fj(Pj)V(j>v) + G(p G ) (31) 

This class of processes includes the decay of the top quark t — > Wb + G in the SM, 
that we will discuss later on. By summing over all final polarizations and averaging 
over initial ones, the square modulus of the amplitude for the process (|3l|) is given, 
in the massless limit for fj, by 

li:\M\^ 3 -^f(M^\ g ,?)F f ( t ,u) (32) 

pol r V 

where the expression for the function Ff(t, u) can be found in the appendix A2. Here, 
the Mandelstam variables are defined as 

t = -^2 (P< ~Pg) 2 - 1, u = fa - p v ) 2 , s = x G -t-u, (33) 
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where p\ = Mj., p v = My, x v = and xq = jjft- By the same procedure 

fi fi 

explained in the previous section, we got the total width r(/j — > fj V + Gx) for the 
inclusive KK graviton production. In the massless limit for fj, this is given by 

(34) 

where GV = g 2 /(A^2M v ), and 

/•(l-v^ 7 ) 2 5/2-1 /'( 1_ v /5r v 7 ) 2 /"*+ 

^/(^a,^v^) — / (^g) / / dtFf(t,u) 

= !i^£ (1 _ M _ av)(1±A) _ 1 + av A= ji_ (35) 

Ztl y (1 — M — XyJ z 



2 



with x A = -j^. 

/ i 

In table 1, we present the numerical results for If{x&,xy,§), at xa = iy = 0, 
and in fig. |2| we plot the function 1/(0, xy, S)/If(0, 0, 5) versus Xy, for representative 
values 6 = 2,3, ... ,6. In fig. |3|, we plot the differential widths defined in eq. (p0|) 
versus = v^^g, and for 5 = 2, 3, . . . , 6. We consider both the massless vector 
boson case and the massive case with ry = 0.45. This value of ry is relevant in 
the top quark decay, where rw = Mw/M t ~ 0.45 (M t is the top quark mass). By 
comparing the distributions of the vector boson and fermion decays in fig. ^, we see 
that in general the positions of the maximum for the fermion distributions are closer 
to zero than in the vector boson case. 

By requiring unitarity for the perturbative expansion, 

m^fiV) <l (36) 

where the total width T(fi — > fj V) (in the fj massless limit) is 

m - fjV) = ('^l 2 + l^l 2 ) PM (37) 

and p(x) = 1 — 3x 2 + 2x 3 , we obtain 

^<Mn( II 32 " M ll f. (38) 

\If(x A ,Xy,d) Ss-lJ 

From the If values in table 1, we get for xy = x A = the following limits from 
unitarity 

M f <M D x (5.5, 4.9, 4.4, 4.0, 3.8) , (39) 
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where the numbers inside parenthesis correspond to 6 = 2, 3, 4, 5, 6, respectively. 

Now we apply these results to the specific case of the top quark decays, where 
V = W ± and gy = —Qa = ^75- The top total width (at tree level, and neglecting 
CKM nondiagonal decays) is, in the b massless limit 



M?G F \V t 
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T(t^Wb) = ^^J^p(x w ) (40) 

where p(xw) = 0.887, and is the standard CKM matrix. Then the total inclusive 
BR for any 6 is given by 

In the case 6 = 2, and xa = 0, we obtain 

BR( ( ^6 +Gx H^ig (42) 

being l f (0,x w , 2) = 8.38 x 10~ 3 (for M t = 175GeV). 

It can be interesting to compare this value, with the rates expected for other 
rare top quark decays both inside and beyond the standard model, also considering 
the potential of future accelerators in this field p2j. In case of negative searches for 



this signal, one will impose an experimental upper bound on the BR of this decay: 
BR(t — > W b + Gx) < where A*°p is related to the experimental sensitivity on 

the top branching ratio. Then, one finds, for 6 = 2, 

M D > (A*p)~* 0.22 M z (43) 

where we expressed the mass scale M t on the right hand side through Mz- Then, we 
can compare this result with the corresponding one for the Z decay, obtained from 
eq.fl27|) for 6 = 2, 

M D >(AfJ" i 0.19M z (44) 

where BR(Z — > // + Gx) < Af xp is assumed. We see that, assuming (at present, 
quite unrealistically) a comparable sensitivity on the two BR's, the lower bounds on 
My) obtained from the gravitational Z and top decays turn out to be comparable, 
too. 



5 Higgs boson decays 

In this section, we analyze the gravitational emission in the Higgs boson decays into 
either two massive gauge bosons or two fermions, 

H(p H ) - V( Pl ) V(p 2 ) + G(p G ) 
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H{jph) - f(pi)f(P2) + G(p G ). (45) 

Using the interaction vertices given in section 2, we obtain the following expressions 
for the square modulus of the amplitudes summed over polarizations 

s^—i 2 = <0k F » (t ' u) 

X 2 M 2 

£l^W/ + d 2 = Nf-L^Fi^u) (46) 

pol 6 1V1 P 

where the functions F^^(t, u) can be found in appendix A2. The definition of Man- 
delstam variables t, u, and s is given by 

t = i (pi + Pg) 2 - ar/,v, w = i (Pa + Pcf ~ Xf,v, s = x G -t-u (47) 

H H 

M 2 

with Xi = jjt, i = V, F. Then, the inclusive total widths are given by 



M„X}S S ^ (M H \ M , 



T(H^SS + G X ) = N,^—^) / / (lA|I/|f) (49) 
where 

/^(za,^/,^) = | dx G (xg)^- 1 I tit y duF^ f (t,u). (50) 

Here, the integration limits are the same as in V — > ff + G (see eq. (|i"8|) ), with 
i/ ^ iy in the case of 1^. The coefficient k in eq.(|4*5|) is equal to 1, unless the two 
final vector bosons are identical particles. In the latter 

The tree level decay widths of H — > ff and H — > VV are given by 

m - / /) = Nf¥j^-p f (x f ), T(H ^VV) = K^^pvixy) (51) 
o7T 8 7rv2 



where pf = (1 — Ax) 2 and pv{x) = y/l — 4x (1 — Ax + 12x 2 ). 
The unitarity conditions here require 

m b <u d ( ™>>><>*>£ f, m h <mJ [V ^m;> \* 

\r H (x A ,Xf, 5) S 8 -i J \Ih{xa,x v , 5) Ss-i) 

for H — > ff + G and H — > VV + G, respectively. In the limit x 2 , Xy — > and ia = 0, 

we get 

M H <M D x (5.5, 5.0, 4.5, 4.2, 3.9) , M H <M D x (5.6, 5.0, 4.6, 4.2, 3.9) , (53) 



15 



for 5 = 2, 3, 4, 5, 6, respectively. 

As an application of our results, we can consider two representative scenarios in the 
SM: the light (M H < 2M W ) and the heavy (M H > 2M t ) Higgs boson. In particular, 
we set Mh = 120 GeV and Mh = 500 GeV, respectively. Then, approximating the 
total H width by the dominant tree-level T(bb) and T(tt+ WW + ZZ), respectively, 
the gravitational decays BR's are given by 

• Light Higgs (Mh < 2M W ) 



S5-1 




2+6 



DR(H-/,/H-(7) = ( i92 — ^ K^) //;(,,.. .■ (51) 



• Heavy Higgs (M H > 2M t ) 

Bm ^WW + G) = [^-J[^LJ I«„(^,x w ,5) (55) 

Ss-i \ ( M H \ 2+S T v i 



( Ss-l \ ( M H \ 2+& T y , 



BR(H^ZZ + G) = IUx a ,xz,S) (56) 



BR(H —fit + G) = (^)(^) Hi**,*,!) (57) 

where the variables Xi are defined cLS Xj — Mf/M H , with i = b,t } W } Z, and A H = 
12x t pf(x t ) + 2p v (x\y) + Pv( x z)- Note that, in eqs. (f48"D and ([51]), Gy ~ > Gp for the 
H decaying both into W's and into Z's. 



In the case 5 = 2, we obtain 



2.2 x 10" 7 
M^(TeV) 



BR(H -*bb + G x ) = ttutk^ > M ^ = 120 GeV ( 5 * 



2 1 x 10~ 5 

BR(F^^^ + G X ) = - , M„ = 500 GeV (59) 

BR(iJ —> Z Z + G x ) = 8 J 4 ^tvV M ^ = 500GeV (60) 

M£(TeV) 

BR(tf^tt + Gx) = ^TeVV M " = 500 GeV ( 61 ) 

where we used, for M H = 120 GeV, I f H (0,x b ,2) = 0.312, and, for M H = 500 GeV, 
i^(0,x w ,2) = 0.0393, 1^(0, x z , 2) = 0.0321, and I f H (0,x t ,2) = 0.00876. 

We can see that, in order to constrain Mjj in the range of a few TeV's for 5 = 2, 
we need a sensitivity on the Higgs BR's of order O(10~ 7 ) and 0(1O -5 ), for the light 
and heavy Higgs boson, respectively. Higher sensitivities are needed to explore the 
case of a larger 5. Such sensitivities on the Higgs BR's are beyond the reach of any 
presently planned experiment by a few orders of magnitude. Anyhow, they could 
become of some interest for physics that might be studied at a Higgs boson factory 
in a not-near future (see, e.g.,p3]). 
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6 Conclusions 



In this paper, we studied the effects of quantum gravity propagating in large ex- 
tra dimensions in a few favoured decay channels of heavy particles. In particular, 
we analyzed the inclusive radiative emission of Kaluza-Klein spin-2 gravitons in the 
following decay channels: the two-fermions decays of massive gauge bosons, heavy 
quarks, Higgs bosons, and the two-massive gauge bosons decay of Higgs bosons. 
Due to the huge number of KK gravitons radiated, the inclusive widths, for a 

/ , \ 2-\-S 

particle of mass M, is only suppressed by a factor of order [jj^J , versus the 

usual factor (j^j arising in quantum gravity in 3+1 dimensions. If the mass of 
the particle is pretty close to the Plank mass in D-dimensions Mp, the quantum 
gravity effects might sizeably affect the heavy particles decays. In scenarios where 
the SM fields propagate in extra dimensions with a ~ TeV compactification scale, 
good candidates for the decaying heavy particles might be the KK excitations of the 
usual SM particles. 

In this framework, we provided analytical results for the square modulus of the 
amplitudes, and numerical results for the inclusive widths. Final-state masses have 
been taken into account, apart from the case of a heavy fermion decay, where a mass- 
less final fermion is assumed. Since, experimentally, the KK gravitons are indirectly 
detected by measuring missing energy and mass in the decay process, we presented 
plots for the distributions of the widths versus the KK graviton mass. We showed that 
the position of the maximum for each distribution is quite sensitive to the number of 
extra dimensions. 

We also discussed the validity of the present perturbative approach for heavier 
masses of the decaying particles. We showed that, when the mass of the decaying 
particles is a few times Md, the radiative widths exceed the corresponding tree- level 
widths, breaking unitarity. While there are not unitarity problems in the effective 
theory for M < Mp, one should keep in mind that, for larger masses, non perturbative 
effects get in general important. 

As an application of our study, we analyzed the decays Z — ► ff+G, W — » ff+G, 
t — > Wb + G, H — > // + G, and H — > WW + G, that can be of interest at present 
and future experiments. In the case of Z decays, the present sensitivity on the 
BR(Z — > ff + E miss ) at LEP1 can push the lower bound on from this decay 
channel not far from the bounds obtained at LEP2 |21[ . Similar bounds are obtained 
from the top quark gravitational decay, assuming (quite unrealistically, at present) 
that some day the experimental sensitivity on its BR will get close to the Z one at 
LEP1. For the Higgs boson decay, we considered the two representative cases of a 
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light (m H = 120 GeV) and heavy (m H = 500 GeV) Higgs. We showed that, in order 
to set lower bounds on M D of order a few TeV's for 5 = 2, a sensitivity ~ 10~ 7 and 
~ 10~ 5 , respectively, on its gravitational BR's is required. The latter sensitivities are 
definitely a few orders of magnitude beyond the reach of the planned experiments for 
Higgs production and study. 
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Appendix Al 



In this appendix we report the Feynman rules for gravitational interactions which are 
relevant for the processes considered in this article. In particular, by means of eqs. (|4]) 
and (H), we obtain 




where we used the convention that particle momenta (indicated inside parenthesis) 
flow along the arrow directions. The expressions of W and X quantities, corre- 
sponding to 3-line and 4-line interaction vertices, respectively, involving 2-vectors 
(W {v \ X^), 2-fermions (W {f \ X^), and 2-Higgs bosons (W^), are given below. 
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Vector 



W /lual3 = ~ Vaf3 fa ■ fa) + Vapfa/xfav ~ ^favfaot 

+ Vfia (Vu/3 fa-fa- faufap) + (rjnaVuP ~ ^V^Vap) } (62) 

X%lp = -^-gM v (rivaVvfS ~ ^VnvVa^ (63) 
Fermion 

W$ = -~{^(fav + fa»)-V^(fa + fa-M fi )} (64) 

x l f J a = ~WWT g ~ 7q7 ^) (9v + 9Als) (65) 



= -^-X f . (66) 
M P u y ' 

• Higgs boson 

= -^{ kl ^ 2u -^(fa-k 2 -M 2 H )} (67) 
where the symbol p stands for p = ^ a p a . 

Appendix A2 

In this appendix we report the expressions for the functions Fy^it^u), Ff(t,u) } and 
F^^(t,u) appearing in eqs.(|I~2"D, (|32"D, and (|^) respectively. Everywhere, the relation 
s = %c — t — u holds. 

• V -> f f + G 



10 
i=l 



T} +) = -^{-16 + 22t-6t 2 + 22u-2tu-6u 2 + lQx f -32tx f 
+ t 2 xj — 32uxj + 2tuxf + u 2 Xf + 2xq (—7 + 2xf) 

- 2x G (22-7t-7u-32x f + 2tx f + 2uxf)} 

Tl~ ] = -^V, (48 - 56 1 + 13 £ 2 - 56m + 26tu + 13m 2 
3 s 2 L 

- 4 (-28 + 13t+ 13m) x g + 52x 2 g ] 
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T 2 (+) = ^(-l + x / )(34-4x G x / -32xj) 

T " = ^(34-4^^-324) 

T 3 (±) = {7f> («~t)} 

T 4 (+) = -^{24 + 36w + 12w 2 + x G (6- 12 a;/) -56x f -68ux f 

- Uu 2 x f + 3u 3 x f + 32x 2 f + 5Qux 2 f -Au 2 x 2 f 

- 2x 2 G (3 + Qu + 2x f -12ux f + 8x 2 f ^j 

+ x G (-12 - Qu + Qu 2 + 48x f + 18ux f - 15u 2 x f 

- -112x 2 f + 16 uxf)} 

T 4 (_) = --^{l2a; G + 4a; G (-9-6w + 4a;/) 

+ x G (-32 + 6M+15w 2 + 16:r / -16w:r / ) 

+ (u-6) (-12-12w-3w 2 + 16a; / + 4wa;/)} 

T 6 (+) = — (l2 + 3a; G - 18a; 2 (-1 + x f ) - 60x f + 80x 2 - 32x 3 , 
3tw L 

+ a; G (27- 750;/ + 32a;j)} 

T 6 H = --^-{36 + 39x G + 8x 2 G -lUx f -56x G x f + 96x 2 f } 

T 7 (+) = ^-{56x^x/ + a: G (88 - 28x f - 53t x f - 53ux f + 32x 2 f ) 

- 2 (-42 + 10t-3t 2 + 10w + 6tw-3w 2 + 64a;/ + llta;/ 

- 7t 2 Xf + lluxf — lltuxf — 7u 2 Xf — 56x 2 + 6tx 2 + 6-ua; 2 )} 

T 7 ( ~ } = |£ {56 a; 2 . + x G (20 - 53 i - 53 u + 32 a;/) 
s 

- 2 (88 + 47t-7t 2 + 47w-lHw-7w 2 -8a;/ + 6ta;/ + 6wa; / )} 
T 8 (+) = ^-{36 + 12t+12a; 2 ; (-1 + x f ) - 12x f -24tx f 

+ 3i 2 a;/ - 14txJ + a; G (9t + 30x / -9ta; / + 16xj)} 

T 8 ( ~ } = ^-{-60-6t + 3t 2 + 12x G - 160 a;/- 14tx/ 

+ a; G (-18 - 9t + 16a;/)} 

rf) = {if> („„«)} 

= -i{8 + 6t + 6w + 22a;/ + tx/ + MX/ + 28a;J + 12a; G (-l + 2x/)} 

T} } = -^(130 + t + u + 24x G + 28x f ) 
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fjV + G 



10 



) = 



=1 



F 1 = ^{(t + u-2x G ) 2 + x v (-12t + llt 2 -9u 2 + 2Ax G -24tx G 

+ Ax 2 G + 2u (-6 + t + 8x G )) - 4 4 (-4 + 8t + 3t 2 + 3w 2 

+ -u (8 + t — 7 xq) — 16 x G — 7 1 x G + 7 4) 

+ 44 (4 + llt+llu-22x G ) -324} 

F 2 = -^(-32-4^ + 34) (-l-xv + 24) 

F 3 = ^(-1-^ + 24) 

F 4 = — (l2-3w 3 -4a; G -24 + 64 + 4w 2 (-l + 3x G ) 

+ u (2 + Qx G - 154) -2x y (-10-32x G -84 + 3a; G 

+ u 2 (-4 + 3x G ) +u (l4+ 12x G -64)) 

- 2 4 (2 + 6m 2 + 12x g -34- m (7 + 3x g )) 

- 44 (1 + 9m-3x g ) -244} 

if t 2 t 

F 5 = _ 2 + 6x G + 54 + 4 + - (2 + x G )-- (6 + 10x G + 34) 
us y 2 2 v 7 

- i v (2 + t 2 (-l + x G )+4x G + 2x G + x G -ta; G (l + 2x G )) 

+ 4 (-6-2t 2 -4:r G + 4 + ^ (9 + x G )) +2xy (5-3t + x G )-A 

F 6 = ^{a + x g -x 2 g -2x v (2 + 4x g -34 + 4) 

+ 4 (~ 12 + 25 x G - 12 4) + 4 ( 20 - 18 x g) - 8 4} 

F 7 = — f-10w 2 + w (34- lit + 31a; G ) 
6s 1 

- 2 (8-22t + 2t 2 + 36x G -llta; G + 144) 

- 2rr y (6-5t + 3t 2 -6 (1 + m + 3m 2 -7x g ) 
+ 44 (16 + 5* + 5u- 22x G ) - 844} 

1 f o 9 * 2 / 3 x G , ,\ 

F 8 = - <^ -3 + t - -x G + -x G - x G + x v [t - -tx G + — (3 + 4x G ) 
H I 24 \22 / 

+ 4( 9 - 2t )-64} 
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F 9 = -{— (-Qu 2 + u (2 + 15x G ) +6 (lO + x G -2x 

- y (3w (-2 + x g ) + (13-4x g ) x G ) - x\ (7 + 2u)-6x y j 
Fio = i{55 + 2t-w + 12x G + xy (ll + 6£ + 6u- 12x G ) + 8x y } (69) 



V V + G 



10 

FX(t,u) = y.hY 



E\ = — {(l6 + t 2 + 20w + w 2 + 2t (10 + w) -4 (7 + t + w) x G + 4x G ) 

x (l-4xy + 12x y )} 

EX = ^{x 2 G (^-6x v + l3x 2 v ^J +x 2 v (4-5u 2 -16x v + 48x 2 v ) 

+ x G xy (3 + 5 u — 12 xy + 56 x y ) } 
^ = («<-*)} 

F 4 y = — {3- 16 i v -M 3 a;v + 52x y -48 x y + x G (-1 + 6 x y ) 

3 x 2 

+ u 2 (3- llxy + 14x y ) + U + u-24x v -6ux v + 32x y ) 

+ it (6 - 23 x v + 48 x y - 12 x y ) + x G (V ^-^ + 5 x v j 

+ u (-y + 41xy-52x y ) +4 (-2 + 11 x y -23 :4 + 6:4))j 

Hf = {Hi (u~t)} 

Hi = — (3- 24x y + 92x y - 176 x y + 96 x v + x 2 G (i + 2x y + 26x y 



tu I v v " v V2 

+ x G (3 -4xy -44x y + 112x y )} 

# 7 y = -{24 + 19w + w 2 + x G (5-22x y ) +t 2 (l-4x y ) 

- 94 x y - 69wx y - 4w 2 x y + 240 x^ + 114 ux 2 y - 24 x y 

+ t (l9 + 2M-69xy-6wxy + 114x y ) 

+ y (-60 -9w + 252x y + 34wx y -424x y + t (-9 + 34x y ))| 

# 8 y = -{6 + x 2 G (l-6xy)-16xy-7t 2 xy + 20x y + 72x y 

+ t (3-6xy + 4x y ) +x G (-8 + 29x y + 8x y +t + 10 av))} 
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Hf = {ff 8 y («<->*)} 

H^o = l9 + t + u-70x v -9tx v -9ux v + 126xy + x G (S + 22x v ) 

• H -> f f + G 

10 

Fl(t,u) = £ff/ 
i=i 

ff/ = ^(4 + t + M-2x G ) 2 

^ = ±{^ + 2s s / + 16sj} 

if/ = {if/ («<-t)} 

if/ = i{-2(-2-« + 2« )(3(l + «)-(4 + «)x / + x (-3 + 2x / ))} 
ff/ = {ff/ («<-t)} 

H( = ^-{l2 + x G (9-8x f ) -48^ + 324} 

if/ = i {72 + 28t + t 2 + 28w + 2tw + w 2 + 4:^-6 (4 + t + w) x f 
+ 4x G (-U-t-u + Axf)} 

H( = i|6 (3 + t) + (16-7t) x/ + a; G (-!5 + y + 8x/)} 

ff/ = {ff/ («<-t)} 

#/ = 39_5t_5 M _ 14^ 
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Figure 2: Values of the integrals Iy (top-left) , Ip (top-right), r R (bottom- left), and l\ 
(bottom-right) versus ry, ryy, rj, and ry respectively, evaluated at xa = 0, and divided by 
their corresponding values at x% = 0. A superscript /( + ) is understood for Iy. In each plot, 
lower curves correspond to higher extra dimensions, for 5 = 2, 3, 4, 5, 6, respectively. Here, 
ri = jrf, Xi = rf, where M is the mass of the decaying particle. 
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Figure 3: Widths distributions R = - sr ^. — -£^- ; versus rc, for xa = 0, and for 6 
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2,3,4,5,6. Here, rc = r i = 37 with i = f,V and M is the mass of the decaying 
particle. Plots relative to the V — > / / + G and fa — ► faV + G decays are shown. 
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